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Abstract 

We study the shock structures in three-states one-dimensional driven- 
diffusive systems with nearest neighbors interactions using a matrix prod- 
uct formalism. We consider the cases in which the stationary probability 
distribution function of the system can be written in terms of superpo- 
sition of product shock measures. We show that only three families of 
three-states systems have this property. In each case the shock performs 
a random walk provided that some constraints are fulfilled. We calculate 
the diffusion coefficient and drift velocity of shock for each family. 

PACS 02.50.Ey, 05.20.-y, 05.70.Fh, 05.70.Ln 

1 Introduction 

During recent years, much attention has been paid to the study of out of equi- 
librium low-dimensional diffusive particles systems [TJ [5]. These systems are 
kept far from equilibrium by maintaining a steady state particle current. This 
results in various critical phenomena including boundary-induced phase transi- 
tions, phase separation and spontaneous symmetry breaking. For instance the 
evolution of shocks in these systems is an indication for the phase separation 
phenomenon. Macroscopically these shocks are special solutions of nonlinear 
hydrodynamic equations for coarse-grained quantities like the density of parti- 
cles. Nevertheless in order to understand the structure of the shocks and also 
the nonlinear nature of them it is necessary to study the microscopic dynamics 
of these systems in details. 

A simple two-states lattice gas model, which shows a variety of interesting criti- 
cal phenomena including shock formations, is the Partially Asymmetric Simple 
Exclusion Process (PASEP) with open boundaries. In the PASEP the particles 
are injected into the system and also extracted from the boundaries while hop in 
the bulk of a discrete lattice [3] . Recent investigations show that for the PASEP 
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with open boundaries (and even on an infinite lattice) a traveling shock with a 
step-like density profile might evolve in the system provided that microscopic 
reaction rates are adjusted appropriately [U [SJ [U [U] • The shock position in 
this case performs a random walk while reflects from the boundaries. It is also 
known that for the PASEP multiple shocks can evolve in the system provided 
that specific constraints are fulfilled. It turns out that the PASEP is not the only 
two-states lattice gas model in which the dynamics of the single shock under the 
Hamiltonian of the system is a random walk. There are also two other systems 
of this kind i.e. the Branching-Coalescing Random Walk (BCRW) and also the 
Asymmetric Kawasaki-Glauber Process (AKGP). The steady state probability 
distribution functions of these three systems are known to be made of superpo- 
sition of product shock measures [5] . 

In the context of one-dimensional out of equilibrium reaction-diffusion systems, 
the Matrix Product Formalism (MPF) has shown to be one of powerful technics 
in order to study of both the steady state and also the dynamical properties 
of these systems j2]. According to this approach the probability distribution 
function of the system at any time (including the steady state) can be written 
in terms of the matrix element (for open boundary conditions) or trace (for 
periodic boundary conditions) of products of some non-commuting operators. 
At large times these operators are time independent while they can be time 
dependent at finite times. For the systems with nearest neighbors interactions 
these operators satisfy an algebra which can have finite or infinite-dimensional 
representations. For the three above mentioned systems, it has been shown that 
the conditions for the existence of a single traveling shock with random walk 
behavior in the system are exactly those for the existence of a two-dimensional 
representation for their quadratic algebras [10]. On the other hand, it is well 
known that the quadratic algebra of the PASEP has also an n-dimensional rep- 
resentation given that some extra constraints on the boundaries and bulk rates 
are held [TTJ [T^] . It turns out that these constraints are exactly those for the 
existence of n — 1 consecutive shocks with random walk dynamics in the system. 
In this case the steady state probability distribution function of the PASEP can 
be expressed in terms of superposition of n — 1 shocks. 

In this paper we consider the most general three-states reaction-diffusion system 
with nearest neighbors interactions and open boundaries. The system is defined 
on a lattice of length L. Throughout this paper, we assign the letters A, B and 
E to these states. The E's are associated with the empty sites or holes in the 
system while A's and -B's are associated with particles of different types. We 
define the total density of particles as the sum of the densities of particles of 
type A and B at each site. We also assume that the density of A particles in a 
given site is proportional to the density of B particles at the same site. In terms 
of the total density of the particles and the density of vacancies, the three-states 
system is basically a two-states system. In present work, we firstly investigate 
the conditions under which the stationary probability distribution function of 
the three-states system can be expressed in terms of superposition of single 
shocks. Our approach will be based on the MPF. From the MPF point of view, 
one can simply try to map the resultant algebra of the three-states system into 
that of three two-states above mentioned models and find the constraints on the 
microscopic reaction rates under which this mapping is possible. This guaran- 
tees that the resultant algebra of the three-states system has a two-dimensional 
representation which in turns indicates that the steady state distribution of the 
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system can be written in terms of superposition of single shocks. In the case 
of mapping the quadratic algebra of the three-states system into the PASEP's 
algebra, the constraints are sufficient for expressing the stationary probability 
distribution of the system in terms of superposition of multiple shocks. 
Secondly we study the dynamics of a single product shock measures in the 
three-states systems. The dynamics of a single product shock measure in two- 
states systems has already been studied and as we mentioned in the PASEP, 
the BCRW and the AKGP the shock performs a random walk given that some 
constraints are satisfied. For our three-states system the conditions under which 
the quadratic algebra of the system has a two-dimensional representation are 
exactly those for the product shock measure to have a random walk dynamics. 
We obtain the diffusion coefficient and the shock drift velocity in each case. 
This paper is organized as follows: In section [5] we briefly review the basic con- 
cepts of the MPF especially for two and three states system. The quadratic 
algebras of the PASEP, the BCRW and also the AKGP are reviewed in section 
[U We also review the constraints under which their quadratic algebras have 
two-dimensional representations. In section f?] we map the quadratic algebra of 
the most general three-states system (in terms of the total density of the par- 
ticles and the vacancies) to the quadratic algebras of these three models and 
find the constraints under which the mapping is possible. The dynamics of the 
shock is discussed in section [5] The results and conclusion are brought in the 
last section. 



2 The Matrix Product Formalism (MPF) 

In this section we review the basic concepts of the MPF introduced in [3] . Let 
us define P{C;t) as the probability distribution function of any configuration C 
of a Markovian interacting particle system at the time t. The time evolution of 
P(C; t) can be written as a Schrodinger equation in imaginary time 

j f P(C;t) = HP(C;t) (1) 

in which H is a stochastic Hamiltonian. The matrix elements of the Hamil- 
tonian are the transition rates between different configurations. For the one- 
dimensional systems defined on a lattice of length L with nearest neighbors 
integrations the Hamiltonian H has the following general form 

L-l 

H=Y^ h k,k+i + hi + h L . (2) 
fc=i 

in which 

h kik +i = Z®^- 1 ) <g> h <g> J®^-*- 1 ) 

hi'= ftW ®1^ L -^ (3) 

h L =1®^-!) (giftW 

For two-states systems X is a 2 x 2 identity matrix and h is a 4 x 4 matrix 
for the bulk interactions. The most general form for the bulk Hamiltonian of a 
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two-states system in the basis (00, 01, 10, 11) is given by 
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In this basis and 1 stand for a vacancy and a particle respectively. Requiring 



the conservation of probability one should have u>a = ^ 
general form the interaction rates are 



In the most 
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Diffusion of particles W32 , ^23 

Coalescence of particles W34, 1^24 

Branching of particles W43 , W42 

Death of particles ^13, ^12 

Birth of particles ^31, ^21 

Pair Annihilation and Creation ^u, luh 

Injection and Extraction of particles at the first site a, 7 

Injection and Extraction of particles at the last site 5, f3. 



In the stationary state we have HP*{C) — 0. Let us define the occupation 
number ti where r, = if the site i is empty and Tj = 1 if it is occupied by a 
particle. According to the MPF the stationary probability distribution P*({ri}) 
is assumed to be of the form 



P*({n}) = ]r(W\ ^(nD + (1 - n)E)\V). 



(5) 



The function Zl in ([5|) is a normalization factor and can be obtained easily using 
the normalization condition. The operators D and E stand for the presence of 
a particle and a vacancy at each site. These operators besides the vectors (W\ 
and \V) should satisfy the following algebra 



(G) 



The operators E and D are auxiliary operators and do not enter in the calcula- 
tion of physical quantities. 

For three-states systems the Hamiltonian H is given by ((2]). In this case X is a 
3x3 identity matrix, ft, is a 9 x 9 matrix for the bulk interactions and hS l > (h^) 
which stands for the particle input and output from the left (right) boundary 
is a 3 x 3 matrix. One should note that by requiring the conservation law for 
the probabilities the matrix h has only 72 independent elements. The boundary 
matrices and hS r > have also 6 independent elements each. Introducing the 
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basis (00,01,02,10,11,12,20,21,22), in which 0, 1 and 2 stand for a hole, a 
particle of type A and a particle of type B respectively, we have 
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for the bulk Hamiltonian and also 
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for the boundaries. As we mentioned the conservation of probabilities in ([7]) 
requires xa = x ji f° r h j — 1; ' ' ' 1 9- At the left boundary the particles of 

kind A and B are injected into the system with the rates a\ and 02 respectively. 
These particles can also leave the system from there with the rates 71 and 72. 
There are also possibilities for changing the particle types at this boundary 
The particles of type A (B) are converted to the particles of type B (A) with 
rate 03 (73). The same processes can also take place at the right boundary for 
the particles of type A and B i.e. the particle injection (with the rates Si and 
#2), the particle extraction (with the rates [3\ and P2) and also the particle type 
conversion (with the rates £3 and /3 3 ). In the three-states case let us dehne some 
new notation. We introduce two occupation numbers, Tj and 6i, for each site 
of the lattice, where Tj = 1 if the site i is occupied by a particle of type A and 
otherwise. Similarly 61 = 1 if the site i is occupied by a particle of type B 
and otherwise. As the particles are assumed to be subjected to an excluded- 
volume interaction, only one of Tj and 9i can be nonzero and the configuration 
of the system C is uniquely defined by the set of occupation numbers {Ti,9i}. 
According to the MPF the stationary probability distribution P*({Tj,0j}) for 
an open system is assumed to be of the form 



1 L 

P*({n, Oi}) - —(w\ WnA + e z B + (l-n- ^)e)|v). 

Zj L 



(10) 



i=l 



The function Zl in (JTOJ) is again a normalization factor and can be obtained 
easily using the normalization condition. The operators E, A and B stand for 
the presence of a vacancy, a particle of kind A and a particle of kind B at 
each site. These operators besides the vectors (W\ and \V) should satisfy the 
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following algebra 




As before the operators E, A and B are auxiliary operators and do not enter 
in the calculation of physical quantities. Using the Hamiltonian of the system 
given by {Z])-([3J ancl CU) the quadratic algebra associated with the most general 
three-states reaction-diffusion model can be obtained. In order to calculate the 
mean values of the physical quantities, such as the mean density of particles at 
each site, one should either find a matrix representation for the algebra or use 
the commutation relation of the algebra directly. 

In the next section we will review the two-states systems in which the Hamil- 
tonian of the system develop a single product shock measure as a random walk 
provided that the reaction rates fulfill some constraints. 



3 Shocks in two-states systems 

In [5] exact traveling wave solutions are obtained for three families of one- 
dimensional two-species reaction-diffusion models with open boundaries. It has 
been shown that for these models the stationary probability distribution func- 
tions can be written as a linear combination of Bernoulli shock measures defined 
on a lattice of length L as 

1 - Pi \ ^ I 1-P2 



w=[ P r ) *{ » ) • °- h - L (12) 

provided that some constraints on the reaction rates are fulfilled. In (fT2|) p\ and 
pi are the densities of particles at the left and the right domains of the shock 
position k respectively. The time evolution of (|12[) generated by the Hamiltonian 
of the above-mentioned models is given by 

£\k) = di\k - 1} + d r \k + 1) - (di + d r )\k) , < k < L (13) 

which is a simple single-particle random walk equation for the position of the 
shock k. The shock positions hops to the left and to the right with the rates di 
and d r respectively. In the following we briefly review the PASEP, the BCRW 
and the AKGP and the conditions under which we have an exact traveling 
shock solution for the model [5] . These models are defined on an integer lattice 
of length L, each site of the lattice i (1 < i < L) is either empty or occupied by at 
most one particle. By applying the standard MPF and defining the operators D 
and E as the operators associated with the existence of particles and vacancies, 
one finds the following results for the above mentioned models: 
For the PASEP the non- vanishing rates in are W32, o., (3, 7 and 5. Using 
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© the quadratic algebra of the PASEP has now the following quadratic form 



CJ23DE - W 32 ED = (lV 2 3 - ^32)(D+E) 

(W / |(aE- 7 D) = (W\ (14) 
{PD-5E)\V) = \V). 

The auxiliary operators E and D are real numbers proportional to unity in this 
case. It is known that the Hamiltonian of the PASEP has an eigenvector with 
zero eigenvalue which can be written in terms of superposition of single product 
shock measures. The non-vanishing rates together with these densities should 
satisfy the following conditions 

P2(l~Pl) _ U> 2 a 
pl (1— (O2) "32 

a(l - pi) - 7P1 = - pi)(w 2 3 - ^32) ( 15 ) 

/3p 2 - 6(1 - p 2 ) = P2(l ~ P2)(^23 - OJ32). 

In the bulk of the lattice the shock position hops to the left with the rate 
di = (uj23 — ^32) Pl p^Zp^ and to the right with the rate d r = (0^23 — ^32) ■ 
The shock velocity and the shock diffusion coefficient will then can be read as 
v s = d r — di and D s — (d r + di)/2. It has also been shown that under the 
conditions (jT5J) the algebra (fT4^) has a two-dimensional representation [TU] . The 
n-dimensional representation of (|14| describes the stationary distribution of the 
PASEP with n—1 consecutive shocks. It has been shown that an n-dimensional 
irreducible representation for any finite n exists provided that the following 
constraint is satisfied by the bulk and the boundary rates 



{ ^l ) l-n = K+{a)7)K+{/3)5) (16) 
LO23 



in which we have defined 



, , —u + v + 1 + J(u - v - l) 2 + Auv 

k+(u,v) = . (17) 

2u 

For the BCRW the non- vanishing parameters in ^ are LO24, u>4St ^34i ^43> ^32, 
U23, a, (3 and 7. The dynamics of a single product shock measure under the 
Hamiltonian of the BCRW is a random walk provided that 

1— Pl _ "24 +"34 
Pl "42 +"43 

1— Pl "23 ( 1 

7 = ^~pT' a + (1 _ Pl)^32 - ^^^43 + P1W34 

while the density of the particles at the right hand side of the shock position 
is zero P2 = 0. The shock position hops to the left and right with the hopping 
rates di = (1 — pi)i^32 + P1W34 and d r — U!ia/pi. The shock drift velocity and 
diffusion coefficient are then v s = d r — di and D s = (d r + di)/2. Using the 
standard MPF the quadratic algebra of the BCRW is obtained to be 

EE - EE = 

W23DE + LU24B 2 - (oj 32 + W42)ED = ED - ED 

-{0J23 + ^43)DE + W34D 2 -I- w 32 ED = DE - DE , , 

W43DE - (w 2 4 + ^ 34 )D 2 4- w 42 ED = DD - DD ^ ' 

(W\(aE - 7D) = (W\E = -(W|D 

0D\V) = E|V) = -D|V). 
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It has been shown that the steady state probability distribution function of the 
BCRW defined by ([5]), can be described by a two-dimensional representation of 
(HU) given that the conditions (TBJ) are fulfilled [TU] . 

For the AKGP, the non- vanishing parameters in are W12, ^13, ^42, ^43, W32, 
a and [3. The particles are allowed only to enter the system from the first site 
with the rate a and leave it from the last site of the lattice with the rate j3. 
In this case we have p\ = 1 and P2 = 0. There are no additional constraints 
on the rates for this model. The dynamics of shock measure generated by the 
Hamiltonian of the system will be a simple random walk on the lattice and the 
shock position hopping rates are di — W13 and d r — The shock drift velocity 
and diffusion coefficient can now be easily calculated. It has also been shown 
that the quadratic algebra of the AKGP given by 

W13DE + W12ED = EE - EE 

-(uj\2 + 0J32 + w 42 )ED = ED - ED 

-(w 13 +w 43 )DE + W32ED = DE-DE 

w 43 DE + U42ED = DD-DD { 1 

<W|aE = (W|E = -_(W\T) 

PD\V) = E\V) = -D\V) 

has a two-dimensional representation. The stationary state of the system can 
be written in terms of superposition of product shock measures (|12[) . 
In the following section, by defining the total density of particles, we map the 
quadratic algebra of our three-states system into (fTi|) . (fT^|) and (120|) . We find 
the conditions under which such mapping is possible. 



4 Mapping of algebras 

In what follows we consider the most general three-states system consists of 
two species of particles and vacancies and look for the constraints under which 
the steady state probability distribution function of the system can be written 
in terms of superposition of product shock measures. It turns out that the 
dynamics of such product shock measure under the Hamiltonian of the system 
is a random walk under the same constraints. In the most general case in a 
three-states system, the particles belong to two different types A and B and the 
vacancies which will be denoted by E. We define the total density of particles as 
the sum of the densities of A and B particles at each lattice site. We assume that 
at any given site of the lattice the mean occupation value of the particles of type 
A is always proportional to the mean occupation value of the particles of type B. 
From the MPF point of view, this means that the operator A, associated with 
the presence of an A particle at a given site, should be proportional to the the 
operator B which is associated with the presence of a B particle i.e. A = 
in which r is a real and positive parameter. This implies that the operators 
A and B commute with each other but not necessarily with E which in turn 
means the stationary probability for occurrence of any configurations of type 
• • • EEE AAABABABBBABEE ■ ■ ■ does not depend on the exact position of 
A's and £?'s in each block surrounded by E's rather to the number of them. We 
should emphasis that such constraint does not result in a trivial two-states model 
even at the microscopic level (see [T3] for instance). Now by defining the total 
density of particles operator D as D :— A + B the associated quadratic algebra 
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of the three-states system (ITT|) can be written in terms of the two operators D 
and E and the auxiliary operators E, A and B. We also assume that A = 
and define D := A + B. Under these assumptions one finds a quadratic algebra 
in term of the four operators D, D, E and E. This algebra can be regarded as 
an associated quadratic algebra of a two-states system with open boundaries. 
Now the question is "what do we get if we try to map this quadratic algebra into 
the quadratic algebras associated with the PASEP, the BCRW and that of the 
AKGP?". Such mapping will indeed impose some constraints on the reaction 
rates of the three-states system. The first result is that the quadratic algebra of 
the three-states system can have two-dimensional representations by imposing 
extra constraints. This means that the stationary distribution of the three- 
states system, when it is described in terms of the total density of particles 
and the density of vacancies, can be written in terms of superposition of single 
product shock measures. Moreover, if the total density of the particles has a 
step-function structure, it will evolve similar to a random walker in continues 
time under the Hamiltonian of the system. In the following, we map the algebra 
of the three-states system in terms of D and E into (TH| . ([!§]) and and 
investigate the outcomes. As we mentioned, this will impose some conditions 
on the reaction rates of our three-states system 1.6- X ij 's introduced in © 
and also the boundary rates in ([5]) and The constraints (TT5|) and (fT5|) in 
the case of mapping into the PASEP 's and the BCRW's algebra will then be 
applied to the total density of particles and the newly defined boundary rates. 
These constraints will guarantee that the dynamics of (TT2|) will be a simple 
random walk and that the steady state of the system can be written in terms 
of superposition of these shocks. 



4.1 Mapping into the PASEP's algebra 

Here we introduce the conditions under which the quadratic algebra of the three- 
states system can be mapped into ([T4|) . One should note that the quadratic 
algebra of the PASEP does not contain any quadratic terms of types E 2 or 
D 2 . Since we have defined A = t^—D and B = -r^— D; therefore, the quadratic 

l+r 1+r ' ' " 

algebra of the three-states system should not contain the quadratic terms of 
types A 2 , B 2 or E 2 . This imposes some constraints on the matrix elements 
of j7]). On the other hand, in order to have a quadratic term of type DE 
one should only have the quadratic terms of types AE and BE. All other 
combinations of A, B and E should also be eliminated in the quadratic algebra 
of the three-sates system. By applying these constraints we have found that 
the most general Hamiltonian of the three-states system (J7|) should have the 
following form in order to be mapped into ()14|) 
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in which xa = ^2 



for i = 2, • • • , 9. All of the boundary rates in ([5]) and © 



are non-zero in this case. We assume that the auxiliary operators D and E are 
proportional to real numbers d and e. In this case one finds that d = — e. By 
defining e := 0^32 — ^23 one can see that the algebra of the three-states system 
can be mapped into (|14[) provided that 



rx 6S + r l x m + x 65 
rx S 6 + r 2 x 8 9 + x$ 5 
r(x 96 + x 98 ) + x 95 
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-- rxss 
r 2 x 99 . 



(22) 



One should also define the bulk parameters as follows 



LO23 



W32 



= ^24 
x 44 - 

= x 42 
X22 - 



- rx 27 
rx i7 - 

- rX43 

rx 2 3 -- 



X37 
X77 - 
-- X73 
X33 - 



r X 34 

;X 7i 

\.X 7 2 
\X32 



(23) 



and the boundary rates as 
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(24) 



The conditions in (|T5|1 should now be applied to the total density of the particles 
and the new boundary parameters defined in (|24j). 

In |13j the authors have introduced a three-states system with a Hamiltonian 
similar to (|2ip and shown that the dynamics of a single product shock measure 
under this Hamiltonian will be a random walk provided that the constraints 
similar to f2"2"|) - ([2"4"|) are fulfilled. It should be mentioned that since (|14p has also 
finite-dimensional representations, multiple shock structures may evolve in this 
system under the condition (fit)]) . Mapping the quadratic algebra of the three- 
states system into the PASEP's algebra not only confirms the results obtained in 
[13] but also prove the possibility of evolving multiple shocks in the same system 
which has not been studied before. One should also note that the models studied 
in [TJ and [T5] are special cases of the system studied here. 



4.2 Mapping into the BCRW's algebra 

In order to map the quadratic algebra of the three-states system into (fT^|) we 
define the total density operator D and rewrite the algebra in terms of the two 
operators D and E as before. Nevertheless, the auxiliary operators are not 
assumed to be real numbers is this case. One should note that the quadratic 
algebra of the BCRW does not contain any quadratic term of type E 2 . This 
means that all of the entries in the first column of (fTTj) should be equal to zero. 
All other combinations of A, B and E are allowed; therefore, the most general 
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Hamiltonian should have the following form 
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(25) 



in which, for the conservation of the probability, defined so that the 

sum of the elements of each column is zero. The reason that the entries of 
the first row in ([2"5|) are taken to be equal to zero comes from the fact that it 
gives a linear combination of operators with positive coefficients equal to zero. 
This cannot be true unless the coefficients (reaction rates on the first row) are 
zero. We now can map the quadratic algebra of the three-states system into 
((T§|) provided that the nonzero parameters of the BCRW's algebra are defined 
as follows 



W23 

U>24 

^32 := 

W34 :— 

LO42 '■= 

W43 := 



: £24 + rx 27 = H X M + TX 37 ) 
: TTF(£25 + r{x 2 6 + £28) + r 2 X 29 ) 

r (i+ r ) (£35 + r(x 36 + x 38 ) + r 2 x 39 ) 

■ X42 + f£43 = f (£72 + ^£73) 

: 1^7(^45 + r{x ie + x i8 ) + r 2 x i9 ) 

r(l+ r ) ( X 75 + r(x 76 + X 78 ) + T 2 X 79 ) 

= (1 + r)(x 52 + rx 53 ) = ^(.£62 + rx 63 ) 
1 ^ L (x 8 2 + rx 83 ) = ^(£92 + rx 93 ) 
(1 + r)(x 5 4 + rx 57 ) = ^(.£64 + rx 67 ) 
^(£84 + rx 8 7) = ^jf(x 9i + rx 97 ). 



(26) 



There are also some extra constraints that should be fulfilled 



£22 - ?"£23 
X44 — rx4 7 

£55 - r(x 5e 



£58) 



r 2 £ 5 g 



= ^{rx 33 - 
= l{rx 77 - 
= ^(-£65 

= ~ (-£85 



£32) 
£74) 

- 7"(£66 

- r(xa8 



£68) - ^£69) 

£86) - r 2 x 89 ) 



(27) 



= ^(-£95 - r{x 96 + x 98 ) + r 2 x 99 ). 
The boundary rates in this case should be defined as follows 



: ai(l + r) = a>i + a 2 
(3i-rP 3 + S 3 = T ^,(l3 1 +r(3 2 ) 
71 + a 3 - r73 = 1^(71 + 7-72) 



a 

P 
1 

5:=5 l =5 2 = 0. 



(28) 



The last constraint in (j28|) indicates that there is no particle injection from the 
right boundary. The conditions given by (|18[) should now be applied to the 
total density of particles and newly defined parameters in (|27)|) and l|28J) . These 
conditions and definitions are enough for writing the steady state probability 
distribution function of the model in terms of two-dimensional representation 
of the algebra. 
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4.3 Mapping into the AKGP's algebra 



In the case of mapping the quadratic algebra of our three-states system into the 
AKGP's algebra, the Hamiltonian of the three-states system ([7]) should have 
the following form 



h = 





Xl2 


XlZ 


Xu 








X17 





\ 





-X22 


X 2 3 























x 3 2 


-2^33 























X42 


x is , 










X47 











X52 


X53 


x 54 








X57 











X 6 2 


X63 










X67 











X72 


X73 


X74 








-X77 











X S 2 


XS3 


xsi 








X87 








V o 


X92 


X93 


X94 








X97 





J 



(29) 



Since ([2H|) does not contain any quadratic terms of types D 2 and E 2 , the 
quadratic algebra of the three-states system should not contain the terms A 2 , 
B , AB, BA and E 2 ; therefore, the first, the fifth, the sixth and the last two 
columns of (f2"9"|) are chosen to be zero; however, all other combinations of A, B 
and E are allowed in this case. In order to map the quadratic algebra of the 
three-states model in terms of the operators D, E, D and E into the AKGP's 
algebra given by (|2"0|) the nonzero parameters in (|29|) should be related to the 
nonzero parameters of the AKGP as follows 



W12 

Wl3 

0J32 



W42 :— 



W43 := 



12 i- rx 13 ) 
rxvr) 



l+r ( X 

j+^(xu 
x i2 + rx A 3 
\{X72 + rx 73 ) 
(1 + r)(x 52 + rx 53 ) 
1 ± r (a; 6 2 + rx 63 ) 
i± r (x 82 + rx 83 ) 
1 T ¥(x g 2 + rx 93 ) 
(1 + r)(x 5i 
1 ± r (a; 6 4 + rx m ) 
1 ± r (x 84 + rx 87 ) 
^(xgi + rx 97 ) 



(30) 



- rx 57 ) 



besides the following constraints 



£22 
2:44 



rx 2 3 

rX47 



X33 
X77 



(31) 



in which xu's for i = 2,3,4 and 7 are defined according to the conservation of 
probabilities. For the boundary rates one should have 03 = rj3. On the other 
hand we should define 



= (1 + r)ai = ot\ + a 2 

= 1 -r0 3 + 8 3 = i^(/3i+r/3 2 ) 

= 71 = 72 = 

= Si = s 2 = 0. 



(32) 



These conditions are enough in order to write the steady state of the three-states 
system in terms of two-dimensional matrices which satisfy (1201) . 
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Pr 



PrB 
PrA 



1 k L 

Figure 1: Sketch of a single shock at sites A: on a lattice of length L. The density 
profile of A particles is always assumed to be proportional to the density profile 
of B particles both on the left hand and the right hand sides of the shock 
position. 



Pi 

PlB 
PlA 




5 Dynamics of shock 

Let us represent a product shock measure with the shock position on the site k 
for our three-states system as follows 



\k) = 



in which piA {PrA} and pis {Prs) are the densities of A and B particles on 
the left (right) hand side of the shock position respectively. According to our 
assumption that the density of A particles is always proportional to that of B 
particles, one should have pis — rpiA and p r s — rp r A- Now defining the total 
density of particles on the sides of the shock position as pi and p r , we have 
pi = (1 + r)piA and p r = (1 + r)p r A- Sketch of such shock measure can be seen 
in figure [TJ In the first case (mapping the quadratic algebra of the three-states 
system into that of the PASEP) the shock (l33|) has a random walk dynamics 
under lj!?Tj) provided that the constraints P^l) - ([M)l are satisfied. In this case we 
obtain 

Pi = T— " — ( V (34) 

- - K+iP > S) (35) 



" 1 + k + (P,6) 

in which we have used (| 1T[) and the boundary rates a, (3, 7 and 5 arc defined 
according to (124|) . These boundary rates besides the total density of particles 
satisfy the following constraint 

Pr(l-Pl) _ U23 
Pr) "32 

a(l - pi) - jpi = pi(l - pi)(u23 - ^32) ( 36 ) 

(3p r - (5(1 - p r ) = p r (l - p 2 )(^23 ~ ^32)- 

The shock position hops to the left and the right with the rates ^~ Pl W32 and 
^wz2 respectively. 

In the second case (mapping the quadratic algebra of the three-states system 
to that of the BCRW) the shock (|33|) performs a random walk again under |25|) 
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provided that ([2"uT) - (l2"5[) arc fulfilled. The total density of the particle on the 
right hand side of the shock is zero. The shock position in this case hops to the 
left and to the right with the rates (1 — pi)u)32 + Pi^a and W43/ ' p\ respectively in 
which LO32, W34 and W43 are defined in (|26|) . The total density of the particles on 
the left hand side of the shock pi, the boundary parameters a and S, besides the 
non-zero parameters defined in (|26p should now satisfy the following constraints 

1 — Pi _ "24 + "34 

Pi "42 +"43 

1 — Pi "23 fi7\ 

Pi ~ "43 _ 

7 = ^pT"" + ~ Pl)^32 - + Pi^34- 

In the third case of mapping the quadratic algebra of the three-states system 
to the algebra of the AKGP the shock ([55)) performs a random walk under the 
Hamiltonian (|29p . In this case the total density of the particles of the left hand 
side of the shock position is equal to unity while it is zero on the right hand 
side of the shock position. The shock position hops to its leftmost (rightmost) 
site with the rate 1013(1043). There is no need to introduce any extra constraints 
more than those in ([3Tjl in this case. 



6 Concluding remarks 

In this paper we have considered the most general Hamiltonian for the Marko- 
vian three-states systems (two species of particles besides the vacancies) with 
open boundaries. From the MPF point of view, the quadratic algebra of this 
system is generated by three operators A, B, and E associated with the ex- 
istence of two species of particles and the vacancies. We have then assumed 
that the density of particles are proportional to each other. By defining the 
total density of the particles as sum of the densities of the particles, this alge- 
bra transforms into a quadratic algebra in terms of the total density operator 
D and E. Now this algebra can be treated as the quadratic algebra of a two- 
states system. We have shown that under certain conditions this algebra can 
be mapped to the quadratic algebras of the PASEP, the BCRW and that of the 
AKGP. Such mapping has two immediate results: The stationary state of our 
three-states system can be written in terms of superposition of product shock 
measures and that the dynamics of a single product shock measure given by the 
Hamiltonian of the system is simply a random walk. 

One can simply check that the conditions which are necessary for mapping the 
algebra of the three-states system into the PASEP 's algebra, are exactly those 
introduced in |13j but obtained from different approach. In [13] the authors 
have introduced a single product shock measure and found the conditions un- 
der which it has a simple random walk dynamics under the Hamiltonian ([14p . 
In this paper we have shown that (|14[) is not the only way one can define the 
Hamiltonian in order to have such property. There are actually two other ways 
which have been introduced and studied here. In the same reference the authors 
have only studied the dynamics of a single shock. We have shown that, at least 
for the case of mapping the algebra of the three-states system into the PASEP 's 
algebra, multiple shocks can evolve in the system. 

One should note that the procedure introduced in this paper is not the only 
way one can write the quadratic algebra of a three-states system in terms of the 
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two operators D and E. It is worth studying the case where the operators A 
and B are not proportional but related to each other. Our approach can also 
be generalized to the systems with more than three states at each lattice site. 
This is under our investigations [IB] . 
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